Effect of Equilibration on Primitive Path Analyses of Entangled Polymers 
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We use recently developed primitive path analysis (PPA) methods to study the effect of equili- 
bration on entanglement density in model polymeric systems. Values of N e for two commonly used 
equilibration methods differ by a factor of two to four even though the methods produce similar 
large-scale chain statistics. We find that local chain stretching in poorly equilibrated samples in- 
creases entanglement density. The evolution of N e with time shows that many entanglements are 
lost through fast processes such as chain retraction as the local stretching relaxes. Quenching a melt 
state into a glass has little effect on 7V e . Equilibration-dependent differences in short-scale structure 
affect the craze extension ratio much less than expected from the differences in PPA values of N e . 

PACS numbers: 83.10.Kn,82.35.Lr,02.70-c,61.41.+e 



I. INTRODUCTION 

Much of the physics of bulk polymeric material is in- 
dependent of chemical detail and arises instead from the 
connectivity of the constituent chain molecules, provided 
the material is sufficiently dense and the chains suffi- 
ciently long. Since chains can not pass through each 
other, they are confined to "tubes" by topological con- 
straints that are often associated with discrete "entan- 
glements" . The tube model of Edwards and de Genncs 
has been a focus of extensive study over the past few 
decades [lj, and qualitatively describes many aspects 
of entangled polymer dynamics. Perhaps its most famous 
prediction is the relation T) oc N 3 for N 3> N e , where r\ 
is the melt viscosity, N is the degree of polymerization, 
and N e is the "entanglement length" . 

A key construct of the tube model is the "primitive 
path", the shortest contour into which a chain molecule 
fixed at its ends may contract without crossing any other 
chains 1]. Everaers and collaborators have recently 
used this construct to perform "primitive path analysis" 
(PPA) of model entangled polymeric systems [3|.|j]. Their 
results for N e agree with the predictions of a chain pack- 
ing model [j| that explains trends in experimental results 
for melts and semidilute solutions of linear homopolymers 
H H II II H The PPA technique is a promising tool 
for obtaining information about the behavior of individ- 
ual entanglements which has not been accessible through 
other theoretical or experimental means. 

In this paper, we address the relation of melt equi- 
libration to the entanglement length measured by PPA 
analysis. Experimental studies of linear homopolymer 
melts have shown that N e is controlled by the Gaussian 
statistics of chains at large scales 0, H, 0- One might 
reasonably assume that values of N e from PPA analysis 
of model melts would also be controlled by large-scale 
chain structure. Our results, however, show that this 
assumption can fail badly. Well-equilibrated states pro- 
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duced using a double-bridging algorithm |icj are com- 
pared to states produced using the common fast-pushoff 
method. For typical equilibration times, use of the fast 
pushoff produces chain stretching at short scales while 
preserving large-scale structure |lf| . This stretching can 
produce values of N e that are two to four times smaller 
than those for well-equilibrated states, showing that fail- 
ure to accurately equilibrate short-scale chain structure 
can cause large systematic errors in PPA analyses. 

Our values of N e for the well-equilibrated states agree 
with values from previous PPA analyses 0, H and 
from other methods [ljj • We verify that the reduction of 
N e for poorly equilibrated states reflects a real excess of 
topological constraints, and show that early stages of the 
evolution of N e towards equilibrium can occur through 
fast topology-changing processes familiar from tube the- 
ories [2j, such as cooperative chain retraction and con- 
straint release. The changes in N e with improving equi- 
libration are correlated to changes in chain structure. 

We also perform PPA analyses of glassy states. In in- 
terpreting the mechanical properties of polymer glasses, 
it is often assumed that the value of N e in the glass is 
inherited from the melt [H While PPA values for 
N e drop slightly upon cooling from a melt to a glass well 
below T g , the changes are comparable to systematic and 
statistical uncertainties. Glass values for N e can not be 
compared to rheological measurements, but can be ob- 
tained from the extension ratio during craze formation 
|T3. ITU . We find that this measure of N e is less sensi- 
tive to equilibration than PPA values, perhaps because 
deformation into the craze structure removes excess en- 
tanglements produced by local stretching. 



II. POLYMER MODEL: METHODS OF MELT 
EQUILIBRATION 

We employ a coarse-grained bead-spring polymer 
model 16, 17] that incorporates key physical features of 
linear homopolymers such as covalent backbone bonds, 
excluded-volume interactions, chain stiffness, and the 
topological restriction that chains may not cross. All 
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monomers have mass m and interact via the truncated 
and shifted Lennard-Jones potential: 

Ulj{t) = 4e(((^) 12 -(f) 12 )-((f) 6 -(f) 6 )), 

(1) 

where r c is the potential cutoff radius and ULj(r) = 
for r > r c . Unless noted, We express all 

quantities in terms of length a, energy e, and time tlj — 
y ma 2 / e. 

Covalent bonds between adjacent monomers on a chain 
arc modeled using the finitely extensible nonlinear clastic 
(FENE) potential 



U. 



kR 2 

niNE(r) = --pin(f-(r/R ) 2 ), 



(2) 



constant Cqq and Kuhn length Ik are related to the chain 
statistics by 



L-oo — — 



l K 1+ < cos{0) > 



l 1- < cos(9) > ' 



(5) 



Lennard-Jones interactions for nonbonded monomers 
cannot be introduced immediately after creating the ini- 
tial chain configurations because chains spatially overlap 
[23. Instead, a soft repulsive potential is used to re- 
duce the chain overlap through gradual introduction of 
excluded volume interactions. The specific soft potential 
used is 



Usoft{r) 



A[l 



COS I 





r < 2 1 / 6 a 
r > 2 1 / e a 



(6) 



with the canonical parameter choices Rg = 1.5a and k = 
30e/a 2 The chains contain N monomers and have 

an equilibrium bond length Iq = .96a. Chain stiffness is 
modeled using the potential 



Ubend(r) = hend 1 



\bi-i\\bi 

= hend(l ~ COs(0i)), 



(3) 



where f\ denotes the position of the ith monomer on a 
chain and 6i is the angle betweeen consecutive bond vec- 
tors h-i =n - fi-x and h = f i+1 - fj. 

Newton's equations of motion are integrated with the 
velocity-Verlet method 0] and timestep St = .005tlj — 
.012t£j. The system is coupled to a heat bath at temper- 
ature T using a Langevin thermostat 19] . Unless noted, 
the simulation cell is a cube whose size is chosen so that 
the monomer number density is p = 0.85a -3 . Periodic 
boundary conditions are applied in all three directions. 
At this density and T = l.Oe/kB the system is a melt well 
above the glass transition [13 • For this model, published 
entanglement lengths vary from about 70 for fully flexi- 
ble chains (kh P „,d = 0) to 20 for semiflexible chains with 
k bend = 2.0e @,H El mill US- For the chain lengths 
used in this study (N = 350 and N = 500), larger values 
of kbend produce nematic order pH . 

We contrast two common methods for creating initial 
states for the PPA analysis. Both begin by generating 
initial configurations of N c h chains without considering 
excluded volume |22| . Each initial chain configuration 
is a random walk of TV — 1 steps with the bond angles 
chosen to give the desired chain statistics 01 



< R 2 {n) >-- 



nil 



1+ < cos(0) > 



I- < cos(9) > 
2 < cos(9) > (I— < cos{d) >") 

n (1- < cos{6) >) 2 



(4) 



where < R 2 (n) > is the average squared distance between 
monomers separated by chemical distance n. Large scale 
measures of chain structure such as the chain stiffness 



for nonbonded beads. The value of A is linearly increased 
from 4e to fOOe over a time 20t£j, as in Ref. [iJJ]. Unfor- 
tunately, this "fast pushoff" procedure creates significant 
distortions in the chain statistics on length scales compa- 
rable to the tube diameter |l0|, |23j . Chains are stretched 
well beyond the form of Eqn. (0} at intermediate chemical 
distances [Toj . 

After the fast pushoff is completed, the two meth- 
ods of equilibration differ. In the pure molecular dy- 
namics (PMD) method, normal Lennard-Jones interac- 
tions are activated, and integration is continued for up 
to 4.5 • 10 5 txj. This maximum equilibration time is still 
much shorter than the longest relaxation time of the sys- 
tems we consider, which is the disentanglement time Td- 
Only after an equilibration run of more than tw d o the 
chain configurations come to full equilibrium but 
equilibration runs much shorter than are commonly 
used in PMD simulations |ToL Il7| . We therefore refer to 
our PMD-prepared states as "poorly equilibrated" . 

The other equilibration method used is the double- 
bridging-MD hybrid (DBH) algorithm described in Ref. 
[Io| . In addition to standard MD equilibration, Monte 
Carlo moves which alter the connectivity of chain subsec- 
tions are periodically performed, allo wing the chain con- 
figurations to relax far more rapidly 24] . Equilibrated 
chain statistics from extremely long MD runs were used 
to obtain target functions for these Monte Carlo moves 
[Tol IT^] . We therefore refer to the DBH-prepared states 
as "well equilibrated" . 

If a glassy state is desired, we increase r c to 1.5a and 
perform a rapid temperature quench at a cooling rate of 
T = — 2 x 10 _3 (e//cs)/T'Lj- The systems are cooled to 
T = .bbe/ks at constant density, and then to T = Ae/ks 
at zero pressure using a Nose- Hoover barostat 18]. The 
resulting glasses have density p ~ 1.02a -3 . Other tem- 
perature and pressure protocols give similar results for 
N e . 

Once preparation of the system is complete, we per- 
form primitive path analyses as in Refs. @,|4[. Except as 
noted in Section IlIII the procedure is nearly identical to 
that of Ref. Q , and we refer to it as the "standard" PPA 
procedure. All chain ends are fixed in space and several 
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changes are made to the interaction potential. Intrachain 
excluded-volume interactions are deactivated, while in- 
terchain excluded-volume interactions are retained. The 
covalent bonds are strengthened by setting k = lOOe, and 
the bond lengths are capped at 1.2a to prevent chains 
from crossing one another [J. Note that we do not at- 
tempt to preserve self-entanglements, since their number 
is negligibly small for the systems considered here Q| . For 
semiflexible chains, the bond-bending potential is deac- 
tivated by setting kf, enc i — 0. The system is then coupled 
to a heat bath at T = .OOle/fce so that thermal fluc- 
tuations are negligible, and the equations of motion are 
integrated until the chains minimize their length. This 
typically requires from 500 to 2000 tlj- Other variants 
of the PPA procedure are discussed in Section lill Al 

Once the chain contour lengths have been minimized, 
we use the formula given in Ref. |3( to calculate the en- 
tanglement length: 



N K = 



< R„„ > 



(N — 1) < b p 



(7) 



where < R 2 ee > is the average squared end-end distance, 
and < bpp > is the mean bond length at the end of the 
PPA run 0- We also calculate the rms variation of b pp 
for each chain, and report the mean of this quantity as 
Ab pp . An alternative method for calculating N e is to fit 
the primitive path chain statistics to Eqn. Q, with Iq 
and < cos(9) > as fitting parameters Q|. The fit value 
of < cos(0) > is then inserted into Eqn. [S] and N e is 
identified with B- In every case this gave values 
of N e consistent with the values from Eqn. J7J, confirm- 
ing that the primitive paths have Gaussian random walk 
statistics, with N e monomers per Kuhn segment Q. 



III. RESULTS 
A. Dependence of N e on Preparation Method 

Table H] shows results from PPA runs for flexible 
{kbend — 0) and semiflexible (fcbend = 0.75e or 1.5e) melt 
states prepared with the PMD and DBH methods. These 
results are for N c h = 500 chains of length N = 500, so 
finite-size effects are small 0. The measured entangle- 
ment lengths depend dramatically on equilibration proce- 
dure. For states prepared using the DBH method, our re- 
sults for N e agree with values from Ref. Q • N e is smaller 
by a factor of two to three for the poorly-equilibrated, 
PMD-prepared states. 

Another difference we find is that Ab pp is larger for 
poorly equilibrated initial states. Fluctuations in b pp in- 
dicate that friction between chains has prevented stress 
equilibration along the chains. This alone would not 
change the value of N e from Eqn. @. However, N e would 
be decreased if friction prevented the chains from mini- 
mizing the total contour length, for example by trapping 
free loops along chains. 



TABLE I: Primitive path analysis results for flexible and semi- 
flexible melts - standard procedure. All of the PMD states 
listed were equilibrated for 480tlj following the fast pushoff. 
The quoted uncertainties in N e are the errors on the means 
of the distributions of Rl e /((N — l)6p P ) for the N c h chains, 
with each chain considered an independent measurement. 



Method 




N e 


<b pp > /a 


Abpp /a 


DBH 





73.2 ± 2.4 


0.1543 


0.0037 


PMD 





24.9 ± 0.9 


0.2695 


0.0144 


DBH 


0.75 


45.1 ± 1.5 


0.2135 


0.0042 


PMD 


0.75 


18.3 ±0.6 


0.3397 


0.0304 


DBH 


1.5 


28.1 ± 1.0 


0.3078 


0.0119 


PMD 


1.5 


13.1 ±0.5 


0.4510 


0.0612 



To test the potential magnitude of such effects we ex- 
plored different algorithms for obtaining the primitive 
paths. These included beginning the PPA procedure 
at T = 1.0e/fcs and cooling slowly, reducing the ex- 
cluded volume interaction between adjacent monomers 
gradually, and increasing the value of a in the LJ poten- 
tial to reduce friction between chains. Ref. Q specifies 
A.b P p < .006a as a criterion for convergence of PPA runs 
to a state with uniform bond tension in the individual 
chains. For the standard PPA procedure, this conver- 
gence failed to occur for four of the six systems in Table 
[J] Table ITTI shows results for these states with our altered 
PPA procedure. While improved stress equilibration re- 
duces Ab p p by a factor of two, the values of N e increase 
by at most five percent. Thus it seems interchain friction 
is not responsible for the large differences in N e . 



TABLE II: Primitive path analysis results for flexible and 
semiflexible melts with procedure modified to reduce friction. 
Initial states and uncertainties are as in Tabic [I] 



Method 


kbend / £ 




< bpp > /a 


Abpp /a 


PMD 





26.1 ±0.9 


0.2638 


0.0072 


PMD 


0.75 


18.9 ±0.6 


0.3346 


0.0143 


DBH 


1.5 


28.8 ± 1.0 


0.3040 


0.0067 


PMD 


1.5 


13.3 ±0.5 


0.4478 


0.0286 



The entanglement length of a bulk polymeric system 
is often related to large-scale measures of chain structure 
such as the packing length p [f|: 



A^ 



P : 



P<Ree>' 



(8) 



which is the volume occupied by a chain divided by its 
mean-square end-end distance. Extensive experiments on 
linear, Gaussian-chain homopolymers have shown that 
N e oc pp 3 d, IE 0- The packing lengths are .68a and 
.65a for the DBH-prepared and PMD-prepared fcf, e „^ = 
states, respectively. Based on the experimentally ob- 
served power-law dependence, one might expect that 
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the N e for the two systems should differ by only about 
10-15%. The actual difference in the N e measured by 
PPA analyses is a factor of three, as noted above. Sim- 
ilarly, the differences in N e between well- and poorly- 
equilibrated states for kbend — 0.75e and kbend = 1.5e are 
far too large to be explained by any difference in their 
large-scale structure. 

The differences in N e may, however, be partially un- 
derstood in terms of differences in short-scale structure. 
Figure n contrasts the chain statistics for the PMD and 
DBH flexible melt states, prior to primitive path analy- 
sis. To characterize chain configurations, we find it useful 
to define the Kuhn-length-like quantity 



of order a, i.e. small compared to the tube diameter, in 
both cases. 



h(n) 



< R 2 {n) > 
nl 



(9) 



This quantity was found to be more sensitive to equili- 
bration than the other quantities considered in Ref. [Hj • 
In Figure ^ h(n) for the DBH-prepared state increases 
nearly monotonically with n, displaying a Gaussian limit 
lk{n) ~ Ik = 1.84a for chemical distances n > 100. The 
PMD states show pronounced non-Gaussian behavior. 
Results for Ik (n) rise too steeply at small n, reach a peak 
at n = n max , then drop to essentially the same large scale 
value as the DBH state as n — > N [13 . The differences 
in the short-scale chain configurations are qualitatively 
consistent with the difference in N e ; the straightening 
of chains at small scales in the PMD state is similar to 
th at p roduced by increasing k be nd, which also lowers N e 

sum. 




500 



FIG. 1: Chain statistics ifc(n) for PMD (dashed curve) and 
DBH (solid curve) melt states of 500 flexible chains of length 
N = 500. The horizontal line, lu(n) = 1.84a, illustrates the 
large-n Gaussian limit of the well-equilibrated state. 

Figure [21 contrasts the chain statistics of the primitive 
paths for the PMD and DBH flexible melt states. The 
peak in the PMD statistics at intermediate range is sup- 
pressed and the final curve is consistent with Gaussian 
statistics. This is why we get the same N e from Eqn. 
@ or a fit to the full curve. The change in the PMD 
statistics during implementation of the PPA appears to 
be more dramatic than that for the DBH state. However, 
the relative displacement of monomers at n > N e is only 




0.02 



500 



FIG. 2: Chain statistics ifc(n) for PMD (dashed curve) and 
DBH (solid curve) melt states of 500 flexible chains of length 
N = 500, after implementation of the PPA. 



B. Evolution of N e with PMD Equilibration Time 

The effect of short-scale structure on N e can be fur- 
ther examined by following both during sample equilibra- 
tion. Figure |21 shows how the melt chain statistics evolve 
with increasing PMD equilibration time for a semiflexible 
(k bend = 0.75e) N ch = 200, N = 350 system. After equi- 
libration runs of order 10 3 txj following the fast pushoff, 
the chain statistics show the short scale stretching typical 
of poorly equilibrated states . After an equilibration 
run of one Rouse time, tr ~ 2 • 10 5 tlj 42], the chain 
statistics are well equilibrated at chemical distances up 
to about ten monomers, but still out of equilibrium at 
chemical distances approaching the chain length. Table 
IIIII shows additional information on the pre-PPA melt 
chain statistics, as well as PPA results for this system 
for PMD equilibration times up to 4.5 • 10 5 tlj ~ 2tr. 
Results for well-equilibrated DBH-prepared states from 
Ref. are shown for comparison. The value of N e rises 
rapidly, doubling by tr/A and tripling by ~ tr. Note 
that the values of Ab pp drop with increasing t eq over 
roughly the same timescale, indicating that friction ef- 
fects diminish. 

The changes in N e and Ik described above are much 
faster than the disentanglement time; tr is only about 
Td/5 for this system Thus the loss of entanglements 
must occur through a mechanism which is faster than 
reptation. Chain retraction is a mechanism studied in 
nonlinear-response tube theories Q wherein the unentan- 
gled ends of a stretched chain contract rapidly inwards 
along its tube, thereby shortening the tube and releas- 
ing entanglement constraints. This process is much faster 
than reptation for N/N e 3> 1 because motion of the chain 
ends arises from entropic tension rather than stochastic 
diffusion, and because the chain center of mass need not 
move 2j. Since the fast pushoff stretches chains, chain 
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TABLE III: Evolution of PPA results for a N ch = 200, N = 350 semiflexible (kbend = 0.75e) system with PMD melt equilibration 
time t eq , following a fast pushoff. The Rouse time of these chains is tr ~ 2 • 10 5 tlj [l2T |. Error bars are as defined in Table [I] 
In the bottom row, n max , Ik, and < 7? 2 e > are limiting values for N c h — * oo, while N e , b PPl and Ab pp are from Ref. Q. 



teq/TLJ 


Tltnax 


lk{n m ax) / a 


10~ 2 < R 2 ee > /a 


N e 


< b pp > /a 


Abpp/a 


5- 10 2 


13 


2.51 


7.20 


11 ± 1 


0.434 


0.058 


3- 10 3 


18 


2.39 


6.89 


12 ± 2 


0.405 


0.048 


5- 10 4 


19 


2.24 


6.79 


22 ± 2 


0.300 


0.016 


1.0 ■ 10 s 


21 


2.16 


6.22 


25 ± 2 


0.265 


0.010 


2.0 ■ 10 5 


21 


1.99 


5.55 


31 ±2 


0.228 


0.005 


3.0 ■ 10 s 


26 


2.04 


5.82 


34 ±3 


0.220 


0.004 


4.5 • 10 5 


29 


2.00 


6.48 


43 ±3 


0.207 


0.003 


Equilibrated 


350 


2.16 


7.24 


45 ±3 


0.21 


< 0.006 




1 10 100 350 

n 



FIG. 3: Evolution of melt chain statistics lk(n) with PMD 
equilibration time for a 200-chain N = 350 system with 
kbend = 0.75e at T = l.Oe/fcs. The dashed and dotted curves 
are after 500txj and 2 • 10 tlj (about one Rouse time), re- 
spectively. The solid curve shows chain statistics of a well- 
equilibrated DBH state, and the thin horizontal line corre- 
sponds to the equilibrium Kuhn length, Ik = 2.16a. 

retraction begins as soon as the Lennard- Jones interac- 
tions are activated. As shown in Figure Eland Table ITTT1 
the chains rapidly contract at scales comparable to the 
end-end distance. < R^ e > drops as much as four stan- 
dard deviations below its equilibrium value at t eq ~ tr, 
showing that this contraction is not an equilibrium fluc- 
tuation. 

As most of the chains are contracting in this way simul- 
taneously, many topological constraints are also released 
far from chain ends. The combination of chain retrac- 
tion and constraint release 0] accounts for the observed 
evolution of the chain statistics at small chemical dis- 
tances. Another indication of the presence of constraint 
release is that the relaxation of the primitive path length 
L pp = (N — 1) < b P p >, which is well fit by the tube- 
theory prediction for chain retraction [lj 

LM = U + ±= + L 2 £ -P(-^W^ (10) 



has a decay time which is clearly below tr; t pp ~ 8 • 
10 4 t lj ~ OAtr. Studies with different initial states and 
kbend — also show that entanglement loss is more rapid 
than reptation, with t pp about half of tr. 

As noted in Section [HI full equilibration of N e should 
occur only after the disentanglement time r^, when the 
chains have vacated their original tubes. It is interest- 
ing that the N e = 43 value obtained for the longest- 
equilibrated PMD state is consistent with values obtained 
for well-equilibrated states when clearly the chain 
configurations remain far from equilibrium at large chem- 
ical distances. This suggests that short-range order plays 
the most important role in determining entanglements. 
Confirming this would require following the evolution of 
N e to t > Td, which is beyond the scope of this study. 

C. Primitive Path Analyses of Glasses 

The six melt states from Tabled were rapidly quenched, 
and standard primitive path analyses were applied to the 
resulting glassy states. In every case, the value of N e in 
the glass was close to that in the corresponding melt, 
even though melt values depended st rong ly on kt> en d and 
equilibration. It is often assumed [l3lll4| that glasses in- 
herit the melt value for N e , and our results are consistent 
with this assumption to within statistical and systematic 
errors (~5%). 

Nevertheless, it is interesting that the values of N e were 
lower in the glass than in the melt in all six cases. This 
reduction could in principle arise from an increase in in- 
terchain friction associated with the higher density of the 
glassy states, but the values of Ab pp in Table HVl rule this 
out. All were remarkably close to the values for the cor- 
responding melt states in Table [I] indicating that friction 
is not responsible for the changes in N e . 

The small drop in N e from melt to glass might be asso- 
ciated with changes in chain statistics. An affine contrac- 
tion of chains would not change N e . However, while the 
density increases by 15% during the quench, the back- 
bone bond length Zrj decreases by only 0.5%. Chains in 
the glass are therefore stretched at short scales (n < 10) 
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TABLE IV: Primitive path analysis results for glasses. The quoted uncertainties in N e (PPA) are as in Table Values of 
A use the PPA values of N e , with X th = ^NJ /l K . Values of A are measured using the same procedure as Ref. |Tsj|. and 



Ni h = \ 2 l K /lo. 


Method 




N e (PPA) 


Abpp/a 


X th 


A 


Nt h 


DBH 





71.2 ± 2.3 


0.0038 


6.5 ±0.2 


5.9 ±0.4 


59 ±8 


PMD 





23.3 ±0.9 


0.0126 


3.6 ±0.2 


5.9 ±0.2 


62 ±8 


DBH 


0.75 


42.9 ± 1.5 


0.0042 


4.8 ±0.2 


5.4 ±0.8 


58 ± 17 


PMD 


0.75 


16.9 ±0.6 


0.0308 


2.8 ±0.2 


5.2 ±0.4 


58 ± 13 


DBH 


1.5 


26.5 ± 1.5 


0.0121 


3.2 ± 0.1 


3.7 ±0.4 


36 ±8 


PMD 


1.5 


12.2 ± 1.0 


0.0619 


2.1 ±0.2 


4.1 ±0.3 


45 ±7 



relative to an affine contraction from the melt. Reduc- 
tion of N e due to this short-scale stretching would be 
consistent with our results for poorly equilibrated melts. 
Additional entanglements may be created as chain ends 
push outward during the quench. 

Another possible explanation could be that the change 
in density couples to chain-thickness effects associated 
with the PPA. For states of constant density, increasing 
the bead diameter during implementation of the PPA 
decreases the calculated value of N e . The values of N e 
in the glassy states average 94% of those in the corre- 
sponding melts, and (p g iass/ p m eit)~ 1/3 is also 0.94. The 
increase in density from melt to glass could correspond to 
an increase in the effective chain thickness, which would 
be consistent with the observed reduction in N e . We 
leave further discussion of this issue to a forthcoming pa- 
per. 

Values of N e in polymer glasses are often inferred 
from measurements of the plateau shear modulus (G *, = 
ApksT /bN e ) in the rubbery regime just above T g 
Below the glass transition, rheological measurements be- 
come impossible, but N e can be measured via the craze 
extension ratio A = po/Pcraze 01 > where po is the 
density of the undeformed glass and p C raze is the den- 
sity of a stable craze. The assumption that entangle- 
ments act like chemical crosslinks leads to the prediction 
A = X th = WNJoJIk , where Ik is the Kuhn length in 
the glass Experimental and simulation values for 

A at T <C T g are consistent with this prediction for val- 
ues of N e obtained from the plateau modulus 0, 0] . 
However, the prediction has not been tested for direct 
measurements of N e in the glassy state. 

Table E] compares values of X th predicted by values of 
N e from PPA analysis of the undeformed glassy states to 
directly measured values of A. The latter were obtained 
by straining uniaxially at constant velocity and measur- 
ing the ratio of the densities of coexisting uncrazed and 
crazed regions as in Ref. [lj| . The values of X th and A are 
consistent for the DBH states, but not for PMD states 
p5| . Indeed, values of A are relatively insensitive to equi- 
libration, which suggests that the craze extension ratio is 
controlled mainly by large-scale chain structure. Chains 
are greatly stretched as they pass from undeformed re- 
gions into the craze, and this may remove the memory of 



short-scale stretching from the fast pushoff. 

IV. DISCUSSION AND CONCLUSIONS 

We have shown that the density of entanglements in a 
model polymeric melt can depend very strongly on equi- 
libration of chain structure at short length scales. Values 
of the entanglement length for poorly equilibrated states 
prepared using non-crossing chain dynamics (PMD) were 
found to be as much as four times lower than values 
for well-equilibrated states prepared using an algorithm 
which changes chain connectivity (DBH). Coil-packing 
models, which focus on large-scale chain structure, fail to 
predict the magnitude of the differences in N e . Instead, 
the low values of N e for poorly equilibrated states (Ta- 
ble P) can be attributed to local chain stretching caused 
by the fast pushoff procedure used to introduce excluded 
volume interactions. 

At the conclusion of the fast pushoff, chain retraction 
ensues. Values of N e can increase rapidly towards appar- 
ent equilibrium over time scales comparable to the Rouse 
time as short-scale chain structure equilibrates (Tablc lTTTl 
Fig- 13- However, large-scale chain structure remains far 
from equilibrium at tr, and N e may continue to evolve 
until the chains have vacated their original tubes after a 
time Td- Since this time scales as (N/N e ) 3 for noncrossing 
chains, preparation methods which either entirely avoid 
producing the local chain stretching (e.g., through use 
of chain prepacking followed by a "slow" pushoff [i"oj|). 
or accelerate equilibration by altering chain connectiv- 
ity |ToL l24j , are far more suitable for the preparation of 
equilibrated states for primitive path analysis. 

Changes in N e upon cooling from melts well above T g 
to glassy states well below T g were small for all equlibra- 
tion protocols and chain stiffnesses examined. We did 
find that values of N e were uniformly lower in the glasses 
than in the melts; this could be attributed either to chain- 
thickness effects in the implementation of the PPA or to 
the stiffness of backbone bonds relative to intermolccular 
bonds. However, the differences were comparable to the 
uncertainties in our measurements, and our results are 
consistent with the common assumption that the value 
of N e in glassy states is inherited from the melt. Finally, 
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we found that values of N e inferred from measurements 
of the craze extension ratio were much less sensitive to 
equilibration than the values of N e obtained from PPA 
analyses of undeformed glassy states. This may indicate 
that deformation into the craze removes any memory of 
the short-scale chain stretching. 
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